A NON-ERGODIC PROBABILISTIC CELLULAR AUTOMATON 
WITH A UNIQUE INVARIANT MEASURE 



PHILIPPE CHASSAING AND JEAN MAIRESSE 



Abstract. We exhibit a Probabilistic Cellular Automaton (PCA) on {0, 1} Z with a 
neighborhood of size 2 which is non-ergodic although it has a unique invariant measure. 
This answers by the negative an old open question on whether uniqueness of the 
invariant measure implies ergodicity for a PCA. 



1. Introduction 

Consider a random process on where S is a finite set, with local interactions 
and a translation invariant dynamic. There are two natural instanciations, one with 
asynchronous updates of the sites of Z d , and one with synchronous updates. In the first 
case, the model is a continuous time Markov process, known as a (finite range) Interacting 
Particle System (IPS). In the second case, the model is a discrete time Markov chain 
known Probabilistic Cellular Automaton (PCA). 

The relevance of IPS in statistical mechanics, as well as in many other contexts, is 
well established. Let us mention a couple of motivations for studying PCA. First, the 
investigation of fault-tolerant computational models was the motivation for the Russian 
school |13}l5]. Second, PCA appear in combinatorial problems related to the enumeration 
of directed animals [7]. Third, in the context of the classification of (deterministic) 
cellular automata (Wolfram's program), robustness to random errors can be used as a 
discriminating criterion [3]. 

For IPS and PCA, the first question is to study the equilibrium behavior. An equi- 
librium is characterized by an invariant measure, that is a probability measure on the 
state space which is left invariant by the dynamic. An invariant measure [i is attractive 
if, for any initial condition, the state of the system converges (weakly) to fi as time goes 
on. 

By a compactness argument, there always exists at least one invariant measure. There- 
fore, there are, a priori, three possible situations: 

(1) several invariant measures; 

(2) a unique invariant measure which is not attractive; 

(3) a unique invariant measure which is attractive. 

In the last case, which corresponds to the nicest possible situation, the model is said 
to be ergodic. Roughly, an ergodic system completely forgets about its initial condition, 
while a non-ergodic one remembers something forever. 
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A classical foundational question is whether the intermediate case exists. In other 
words, does uniqueness of the invariant measure imply convergence to it ? For monotone 
systems, the intermediate case does not exist. But in general, the question is open. 

For IPS, this question is Open Problem 4 in Chapter 1 of the classical textbook by 
Liggett [8]. In |10j . Mountford proves that the intermediate case does not exist for 1- 
dimensional IPS (that is d = 1). Quoting [10] , "it seems more than plausible that the 
conclusion (...) is true in higher dimensions". However, the question remains unsettled. 
For PCA, the same question is Unsolved problem 3.4-3 in Toom or Unsolved problem 
5.7 in Toom p2]. 

In the present paper, we answer the question for PCA by exhibiting a 1-dimensional 
PCA, model A, corresponding to the intermediate case (Theorem l3.ip . There is a unique 
invariant measure of the form (fiQ + Hi)/2 and the PCA maps ^iq to [i\ and /ii to Re- 
starting from an initial measure fj-o, the probability measure of the state of the system 
is fiQ at even times and \x\ at odd times. Therefore there is no convergence. 

Observe that the situations for IPS and PCA are different: in 1-d, the intermediate 
case exists for PCA, and not for IPS. This is consistent with the situation for Markov 
processes on a finite state space: in discrete time, periodic phenomena may occur which 
result in the existence of the intermediate case; in continuous time, the intermediate case 
does not exist. 

To prove the result for model A, we introduce two auxiliary PCA. The first one, 
model B, corresponds to independently moving particles annihilating when they meet 
(p + p — > 0). The second one, model C, corresponds to independently moving particles 
merging when they meet (p + p — > p). We compute exactly the evolution of the one- 
dimensional marginals for model C (Theorem I5.2|) and models A and B (Prop. 16. 1|) 
starting from a "full" configuration. In particular, it proves that the speed of convergence 
to the invariant measure is of order for the three models. 

Continuous-time versions of models B and C have been studied in the IPS literature 
under the names of annihilating random walks and coalescing random walks, respectively, 
see [UG2E]. The PCA and IPS versions of B and C share the same features: ergodic- 
ity with the invariant measure being the "all empty" Dirac measure, and with similar 
and subexponential speed of convergence. In the IPS setting, the asymptotic speed of 
convergence was given by Bramson & Griffeath [3] for model C, and by Arratia pQ for 
model B. Also, the coupling between the models B and C, that we use in Section 
already appears in Griffeath [61 Ch. 3, Sec. 5] and in Arratia [1] in the continuous-time 
setting. The novelty is that we get exact computations for the PCA models, as opposed 
to asymptotic results for the IPS ones. At last, let us mention that IPS versions of 
models B and C on a finite set of sites have also been studied, see for instance [2] for B, 
[H] for C, and the references therein. 

2. Probabilistic Cellular Automaton 

Let S be a finite set. Denote by A4(S) the set of probability measures on X. Let 
us equip X = S z with the product topology. Denote by M.(X) the set of probability 
measures on X for the Borelian cr-algebra. Weak convergence of (fi n )n to is denoted 
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by n n — > [i. Let K be a finite subset of Z and consider x E £ . The cylinder denned 
by x is the set 

*x* = ju 6 S z , V/c £ K,Uk = ■ 

Given k £ Z and V = , . . . , v n ) G Z", we use the notation k + V for (A; + «i, . . . , k + 
v n ), and the notation V{K) for {i \ 3k £ K, 3v £ V, i = k + v}. 

Let us introduce probabilistic cellular automata, restricting ourselves to 1-dimensional 
models. 

Definition 2.1. The alphabet is a finite set S; the set of sites is Z. The set of con- 
figurations is X = £ z . Given V £ IT 1 , a transition function of neighborhood V is a 
function f : YX — > A4(T>). The probabilistic cellular automaton (PCA) F of transition 
function f is the application Ai{X) — > M(X), fi t-t fiF defined on cylinders by: \/K, 
Vy £ Z K , 

fj,F(*y*) = n f(( x i)iek+v)(yk) ■ 

Let us look at how F acts on a Dirac measure 5 X . The value of all the sites are 
updated. The value x^ of the k-th site is changed into the letter a £ £ with probability 
/((xj)j e fc + y)(a), independently of the evolution of the other sites. 

By specializing Definition 12. T\ we recover two famous models: 

• Assume that V = {0}, then all the sites behave independently. The restriction 
of the PCA to one site is a Markov chain evolving on E. Conversely, any Markov 
chain on a finite state space E can be realized as (a restriction of) a PCA on the 
alphabet E with neighborhood V = {0}. 

• Assume that the transition function / is such that: Vn £ S , f(u) is a Dirac 
probability measure. Then we may view / as a function TX — > EJ. We obtain a 
(deterministic) cellular automaton. 

A PCA F may be viewed as a Markov chain on the state space S z . Thus we borrow 
the classical terminology of Markov chains. 

Definition 2.2. An invariant (probability) measure of F is a probability measure jj, £ 
M(X) such that fiF = u. The PCA F is ergodic if it has a unique invariant measure 
which is attractive, i.e. if 



(1) (i) [3\fi£M{X), nF = n], (ii) £ M(X), vF 




Consider for a moment a Markov chain on a finite state space with transition matrix 
P. Let G(P) be the graph of the matrix P. Classically, we have 

(2) (i) <?=^ Q (P) has a unique terminal component 

(i)+(ii) <^=^ Q{P) has a unique terminal component which is aperiodic . 

In particular, uniqueness of the invariant measure does not imply ergodicity. The sim- 
plest example of a non-ergodic Markov chain with a unique invariant measure is the 
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following: the state space is X = {0, 1} and the transition matrix is 

The unique invariant measure is /i = (5q + Si)/ 2 and for v = Sq, we do not have 
vP n — > /U. 

For PC A, it was an open question to know if (i) implies (ii) in (pQ). The purpose of the 
present paper is to settle the question by proposing a non-ergodic PCA with a unique 
invariant measure. 

To get a hint of the difficulty, consider for instance a PCA F with neighborhood 
V = {0}. Recall that each site behaves independently and as a finite Markov chain P. 
As recalled in ([2]), P may satisfy either [->(«)], [(i),^(ii)], or [(i), (ii)]. We show in Table 
Q]how this gets reflected on the PCA F. 



Markov chain P 


PCA F 


-n(t) 
(«).-"(«) 


-.(0 


(i)i(M') 


(i),(ii) 



Table 1. Finite Markov chain versus "neighborhood PCA". 



Let us justify the Table. If \i is an invariant measure of P, then the product measure 
fj,® 1 ' is an invariant measure of F. Therefore, if P has several invariant measures, the 
same holds for F. Assume now that P is ergodic with unique invariant measure [i. 
One proves easily that is attractive, so F is ergodic. Let us concentrate now on 
the intermediate case for P. If P satisfies then Q(P) has a unique terminal 

component which is periodic, say of period 2. Let (fio + Hi)/2 be the unique invariant 
measure of P. Then F has an infinite number of invariant measures. Indeed, consider 
any (uj)igz G {0, 1} Z , and let (vi)i e z be defined by «j = 1 — Ui for all i. Then the 
probability measure ((di^z^Ui + ®i&iHvi) /2 is clearly an invariant measure of F. 

3. Statement of the main result 

3.1. Model A. Consider the PCA Fa on the alphabet S = {0, 1}, with neighborhood 
V = {—1,0}, and transition function a defined by: 

o(00)(l) = 1/2, o(01)(l)=0, a(10)(l) = l, a(ll)(l) = 1/2 . 

A realization of the Markov chain is obtained as follows. Consider the function 

(4) A: {0,l} Z xW Z ^{0,l} Z 

(xi)iez, (ui)iez ^ (%i)iez , 
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Figure 1 



The transition function of the PCA Fa- 




called the update process. Set U n = (J7i, n )ieZ- Given a {0, l} z -valued r.v. Xq = (Xj ! o)igz ! 
such that U _L Xq, define the sequence of {0, l} z -valued r.v.'s (X n ) n& ^ as follows: 



Then (X n ) ne ^ is a realization of model A. The process U is used to randomly update 
the value of a site, when needed, with — > being interpreted as "keep" and f as "switch", 
and Xi^ n is the state of site i at time n, so that X n = (Aj !n )j S z denotes the state of the 
system at time n. 

3.2. Invariant measure. Let x = (01) z be the configuration defined by: Vn 6 Z, 
X2n = 0, X2n+i = 1- The configuration (10) z is defined similarly. 

Theorem 3.1. The PCA Fa has a unique invariant measure which is fi = (<5(oi)z + 
<5( 10 )z)/2. The PCA is non-ergodic. 

On configurations without 00 and 11, the PCA acts as the translation shift. Therefore 
H = (£(oi)z + ^(io) z )/2 is an invariant measure. Assume that it is the unique one. Then 
fi is non-attractive, the situation being the same as for ([3D: consider v = <5( i)z, then 
vF^ = 5(oi)z if n is even, and vF^ = <5( 10 )z if n is odd. 

The purpose of Sections H] and [5] is to prove Theorem 13.11 



We now define two new PCA, that we call respectively model B and model C . For 
both models, the alphabet is £ = {o, •} and the set of sites is Z. Given a configuration 
u G {o,»} z , the following interpretation holds: if Uj = o, the site i is "empty"; if U{ = •, 
the site i contains a "particle". At a given time step, a particle decides (independently 
of the others and independently of the past) to remain at its site with probability 1/2, 
or to jump to the site on the right with probability 1/2. In model B, if two particles 
collide, then they annihilate. In model C, if two particles collide, they are merged into 
one particle. Let us define the models more formally. 



(5) 



A n+ i — A(X n , U n ). 



4. Two auxiliary models 
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4.1. Model B. It is the Markov chain on {o,»} z defined as follows. Consider the 
function 

(6) B: {o,.} z xW z ^{o,.} z 

with U = {|, ->■}, and 

/• if € {(•o,^U),(o;U t), (••,tt), (••> -►-►)} 

I o otherwise . 

Let ?7 be an update process, defined as in Section I3TT1 Given a {o, •} z -valued r.v. Yq, 
such that U _L lo> define the sequence of {o, »} z -valued r.v.'s (y n )neN as follows: 

(7) y„+i =B(Y n ,U n ). 
Then (Y n ) n ^ is a realization of model B. 



Step n+1 



o — o — o 
• — o — • 
• — • — o- 



proba 1/4 
proba 1/4 
proba 1/4 
proba 1/4 



Step n 



-o- 



Figure 2. The transition function of model B. 



Remarks. In the above presentation, model B is a Markov chain with synchronous 
updates and local interactions, but not stricto sensu a PCA. Indeed, if Yq is deterministic, 
then the r.v.'s Y^i and are not independent, since they are updated using the 

non-disjoint r.v.'s {C/j_i,o, C^,o} an d {Uifi, f7j+i,o}- However, it is possible to give a PCA 
presentation of model Bona larger alphabet. Define the sequence of f{o, •} xU^j -valued 
r.v's (T„)neN by Y n = (Y n , U n ). We have: 

(Y n+1 )i = (B(Y n )i, Ui ,71+1) • 

Thus (Y n ) n is a realization of a PCA on the alphabet {o, •} x U , with neighborhood 
V = { — 1, 0}. The same remark holds for model C below. 

The continuous time version of model B, with exponential holding times, is called an 
annihilating random walk (cf. [61 Ch. 3, Sec. 5]). 
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Step n+1 



-O 



-o — • 



proba 1/4 

q proba 1/4 

^ q ^ proba 1/4 

f q proba 1/4 



Step n 



• — o 



Figure 3. The transition function of model C. 



4.2. Model C. It is the Markov chain on {o,»} z defined as follows. Consider the 
function 

C: {o,.} z xU z -> {o,.} z 

with 

~. _ /* if 0i-l2;,ii;-iUi) 6 {(•o,-»-i/) ) (o«,W f), (••,tt), (••, W)} 

I o otherwise . 

Let U be an update process. Given a {o, »} z -valued r.v. Zq, such that C/ _L Zq, define 
the sequence of {o, •} z -valued r.v's (Z n ) n€ ^ as follows: 

Zn+l = C-{Zm U n ) . 

Then (Z n ) n ^fq is a realization of model C. 

Again, the continuous time version of model C, with exponential holding times, is 
called a coalescing random walk (cf. [61 Ch. 2, Sec. 9]). 

4.3. Links between models A, B, and C. One-step transition of the model B, resp. 
C, defines the mapping 

F B : M({o,.} z )^M({o,.} z ) 

fi i — ► (j,F B , 

F C : M({o,.} z )^M({o,.} z ) 

H i — > iiF c ■ 

<p: { ,l} z ^{o,.} z 



respectively, 



Define 
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with 

= [• if XjXj+i G {00, 11} 
Vi ~\o if Xi x i+1 € {01, 10} 

By extension, define </? : .M({0, 1} Z ) -)■ .A/f ({o, «} z ). 

Lemma 4.1. T/ie diagram below is commutative: 

M{{0,1} Z ) M({0,1} Z ) 



M({o,.f)-^M({o,.f) 

/f (X n ) ng N is a realization of the Markov chain A, then (</?(^n))neN is a realization of 
the Markov chain B. 

Proof. Recall that A and B are defined in (HJ) and ([6|) respectively. We are going to 
prove that: 

(8) <p o A = B o (<p,Id) . 

The statement of the lemma follows. Set 

(xi)i, (ui)i i-A- (xi)i i-^- (ui)i A {m)i ■ 

To obtain ([8]), it is enough to check that yo = yo- This is done by systematic inspection 
in Table [2j Each one of the 32 cases mimicks the commutative diagram: in the first line, 
from left to right, (x_2, X-i, xo)> (w-i,«o)> anc ^ O^-li^o); i n the second line, from left to 
right, (y_i,?/o), («-i,tto), and y = 2/o- 

If the process X is defined by ([5]), relation ([8]) entails that the process Y, defined by 
Y n = ip(X n ), satisfies relation (|7|). □ 

Lemma 4.2. Model B is dominated by model C: for x,u G {o,»} z x U 1, , 

B(x, u) < C(x, u) , 

where < is the coordinate-wise product ordering on {o,»} z ; with o < •. 

Proof. This can be checked directly on the definitions of B and C. Intuitively, particles 
are merged in model C, and annihilate in model B. □ 

Lemma 4.3. The following implications hold: 

[C is ergodic with invariant measure 5 z] 

==^ [B is ergodic with invariant measure 5 z] 

<J=^ [A is non-ergodic with invariant measure (5(oi)z +<5(io) z )/2] • 

Proof. This is a direct consequence of Lemmas 14.11 and 14.21 □ 

Therefore, in order to prove Theorem 13.11 it is sufficient to prove that model C is 
ergodic with invariant measure <5 z. This is the purpose of next section. 
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Table 2. The 32 possible cases. 



5. Model C is ergodic 
Lemma 5.1. Model C is monotone, that is: for z £ {o,»} z ; z £ {o,»} z , u £ , 

z < z C(z,u) < C(z,u) , 

where < is the coordinate-wise product ordering. 

Proof. It can be checked directly on the definition of C. □ 

With this monotonicity, to get the ergodicity, it is enough to prove that <5.zF^ — y 5 z. 
Indeed, consider two realizations of model C, one, say Z = (Z n ) n , that starts with all 
sites occupied, the other, say Z = (Z n ) n , that starts with an arbitrary initial condition, 
their evolution using the same update process U. According to Lemma [57TT at any time 
n £ N, Z n < Z n . 

From now on, we focus on the process Z. Recall that for each n, Z n = (Z^^kez is the 
state of the system at time n. The process Z n is stationary, i.e. invariant by translation, 
since Zq, U, and C are invariant too. Define 

(9) d n = P(Z fc , n = «)=P(Z ,„ = «) . 

This is the density of particles at time n. The density d n can also be viewed as an evalu- 
ation of the distance between Z n and 5 Q z. Indeed, for any finite subset E of Z, consider 
the Hamming distance on {o,«}- E , and denote by Wh the corresponding Wasserstein 
distance on ^({o,*}^). Setting Z E , n = (Zk,n)keE, we have: Wh (Ze,u, S e) = \E\d n . 
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Theorem 5.2. Let T be the time that a simple symmetric random walk on Z needs to 
reach 2, starting from 0. We have 

(10) d n = F(T> 2n) 

r 2n + V 



(11) = 4" n 

n 

In particular, d n ~ 2/y/Wn, hence converges to as n grows. 

In continuous time, when the particles perform a simple symmetric random walk, 
Bramson & Griffeath |3j obtain the same asymptotic behavior for d n , up to a scaling 
factor, as expected. 

Corollary 5.3. Model C is ergodic with unique invariant measure 6 D z. 

Proof. We first prove ([TTj). assuming (fTOjl , Let S = (Sk)k<=N be a realization of the simple 
symmetric random walk on Z, starting from 0. Define M k = maxjS'jjO < i < k}, the 
maximum of the random walk at time k. Recall that T = inf{i > | Si = 2}. We have 

P (T > 2n) = P (M 2n < 1) = 1 - P (Af 2n > 2) 
= l-^P(M 2n >2, S 2n = l) 



1 - P (S 2n > 2) - J> (M 2n > 2, S 2n = £) 



£<i 

According to the reflection principle, for I < 1, P (M 2n > 2, S" 2n = £) = P (5 2n = 4 - £). 
Therefore, 

P(T>2n) = 1 -P(S 2n > 2) -F(S 2n > 3) 
= l-P(5 2n <-2)-P(5 2n >3) 
= P(5 2n G{0,2}) 

--(ax;;,))— r; 1 

Using Stirling's formula, we get 

4 -n/2n + l 



n / \/Trn 



Now let us prove ()10p . Recall that Z = (Z n ) ng ^ is a realization of model C with 
Zo = » z . One can extend the definition of Z via coupling from the past. Consider the 
i.i.d. r.v.'s {U k , n ) (k , n)&xZ with ¥(Uij =t) =V(Uij =-►) = 1/2. For each s £ Z, define 

Z« = (Z< s) ) n > s by 

Z« = . z , Vn > s, Z { n % = C{Z^\{U Kn ) k& ) . 

The starting time of Z^ is s, but, besides that, the dynamic is the same as that of 
process Z. Observe that Z^ = Z. More generally, Z^ s ' has the same distribution as 
{Z- s+n ) n > s . Thus, we have 

d n = P [ Zq q = 
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Figure 4. a. The updating ([4,„)(fc,„) e j- 9 ,o]xM,o]- b - The process 
Z^ -3 ) during the span [—3,0]. c. The process Z^ -5 ) during the span 
[—5,0]. d. The merging of particles. 



In Figure [U we have represented space-time diagrams for the model. The point of 
coordinate (k, n) corresponds to site k at step n. We have also represented the updating 
variables with the following convention: at the point (k,n), there is an arrow pointing 
north if Uk^ n =t and an arrow pointing north-east if U^^n = ~ This allows to visualize 
the evolution of particles in the processes ZW. In FiguresHlb and[Hc, the processes 
and Z^ -5 ) are represented; the grey nodes are the ones whose color depend on updating 
variables outside of the represented window. In Figure HJd, the particles painted in 
orange (gray) are those that merged into the particle present at time and site 0. 

Let I n , n > 0, be the set of indices of particles present at time — n in Z(~ n ) that merge 
into particle at time 0. Either I n = [(i n ,y, a n < b n , in which case Zq Q n ^ = •, or 

I n = 0, in which case Zq = o. We focus on \I n \- For instance, in FigureHld, we have 
G-fnDneM = (1,2,1,2,3,4,3). Observe that 

(12) d n = P (zj Q n) = .) =P(|/ n |>l) . 



We have ao = bo = 0, and |/o| = 1- Define 

L 1. 1, 
P = + + ■ 

We check the following: 

• Assume that a n < b n . Then, 

\a n -l if E/on-i ,-n-i =-> , \b n -l if U bn _ n _i =-?> 

[a n if ,-nr-l =t [bn if t4„ ,-n-X =t 

Thus, |/ n +i| — | in | £ {0, ±1}, and the conditional law of |/ n +i| — \In\ is />• 
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• Assume that a n = b n . Then 



'n+1 



For instance, the third case appears between lines —3 and —2 in Figure Hid. Here 
again, the conditional distribution of |/n+i| — \In\ is P- 
• If I n = 0, then I n +i = 0. 

Consequently (|/ n |)nGN is a random walk with step p, starting from 1, and killed when 
it reaches 0. Using (fT2|) . we obtain (fTUj) . □ 



'0 


if (Ua n 


-1,-n-l 




_ n _l) = (t,->) 


\a n -l,a n - 1] 


if (Ua n 


-1,-n-l 


^„ 


_ n _l) = 


\a n - l,a n j 


if (U an 


-1,-n-l 


^O n 


-n-l) = Kt) 


Ja n ,a n } 


if (Ua n 


-1,-n-l 




_„_!) = (!,!)• 



6. Speed of convergence for models A and B 

Let (A n ) nG N be a realization of model A, with ~ P S -M({0, 1} Z ). The possible 
limits for weakly-converging subsequences of (A n ) ne ^ are of the formp<5( i)z + (l— p)<5(io)z 
for p G [0, 1]. An evaluation of the distance to the limits is given by 

F(A 0}n A hn G {00,11}) . 

Since model A is not monotone, we do not know for which initial measure p this distance 
will be maximized. Hence we evaluate the "speed of convergence" for model A by the 
quantity: 

d£ = max P (A 0>n A 1)n G {00, 11}) . 

«6A1({0,1} Z ) 

The quantity is also the speed of convergence to 5 Q z for model B. Indeed we have 
F(A 0t nA lt n G {00, 11}) =¥(ip(A n ) = •), which implies that 

d£ = max P (5 ,n = •) , 

u£M{{o,.}%) 

where (B n ) n denotes a realization of model B and v denotes its initial distribution 
(B ~ v). 

Recall that d n = 4 -n ( 2ra n l ~ 1 ) is the speed of convergence for model C, see © and 
Theorem 1^21 



Proposition 6.1. We have 



<d^<d n . 



Let (^4 n )n be a realization of model A, with Aq ~ p. If p is the uniform distribu- 
tion on {0, i.e. the r.v.'s A ifi are i.i.d. with P (A 0fi = 0) = P (A 0fi = 1) = 1/2, 
then we shall see that P (Aq^Ai^ G {00,11}) = d n /2. If p = 5-yz, then we have 
P(^4o,n^4i,n £ {00, 11}) = (f n _i/2, which is larger than d n /2. The results can be trans- 
lated to model B: the density of particles at step n is d n /2 if the initial distribution is 
uniform, and it is d n -\/2 if the initial distribution is <5.z. 

The end of the section is devoted to the proof of Prop. 16. 1| through the study of 
model A with the two initial distribution mentioned previously. 
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Figure 5. A realization of process D (gray for g(reen), black for b(lue)). 



We define a new PCA, called model D, which is a coupling of models B and C. The 
alphabet is {o,b, g} and the set of sites is Z. Given a configuration u S {o,6,g} z , the 
interpretation is as follows: if Ui = o then site i is empty; if -Uj = b then site i contains a 
Wue particle; if Ui = g then site i contains a green particle. Particles move as in models 
B and C. When two particles collide, they get merged into one particle as in model 
C. In absence of collision, particles keep their color. In case of a collision, the merged 
particle is colored according to the rules: 

(13) b + b-^g, g + g-^g, b + g-^b, g + b-^b. 

We have represented a realization of model D on Figure [5j The "question mark" nodes 
are the ones whose color depend on updating variables outside of the represented window. 
Define 



■kb -{°,b,g} 
ttc --{o,b,g} 



TT B (o) 
TT C (o) 



o, ir B (b) 

O, TT C (b) 



*,K B (g) 
*,nc{g) 



We keep the same notations for the product applications: ttb ■ {°, b, g} A 
{■n B {ui))i, and ir c ■ {o,b,g} z -> {o,»} z ,(u;); h-> (7rc(«i))i- 



Lemma 6.2. // (D n ) n is a realization of model D, then (7r_B(-Dn))n is a realization of 
model B, and (irc(D n )) n is a realization of model C . As consequences, model D is 
ergodic with unique invariant measure 5 Q z, and d^ < d n . 

Let (D n ) n be a realization of model D with Do being defined as follows: the r.v.'s 
Difl are i.i.d. with P (A^o = b) = P (Z?o,o = 9) = 1/2- At step n, the colors of the 
remaining particles are still i.i.d. and uniform : whatever the shape of the binary tree 
of coalescences leading to the presence of a particle at a given position at time n (see 
Figure [5] for an example) , if the colors of the initial particles are independent and if one 
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of these particle's color is uniform, then, due to (|13p . the resulting color will still be 
uniformly distributed. Therefore we have 

P (7TB(A),n) = •) = P (A),n = b) = ^P (Z> ,„ G {6, <?}) = ip (7T (A),„) = •) = T^n , 

where the last equality follows from Theorem 15, 21 

Now let (D n ) n be a realization of model D with Dq = b z . Define E = (Ei)i by 

E=\ h [i ^> 1 = b 

\g if A,i = o or g 

The r.v.'s are i.i.d. with P (E = b) = F (E = g) = 1/2. Let us justify this 

point. The state at time 1 of a realization of model A that starts from O z is uniformly 
distributed by definition. Hence, the state at time 1 of a realization of model B that 
starts from » z is uniformly distributed. And E has the same law as the latter up to the 
transformation if>i, g -R- o. 

So we have E ~ Do- Observe also that itb(Di) = ttb(E). We deduce that, for all 
n > 1, we have nB(D n ) ~ 7r#(.D n _i). In particular, using flT3J), 

P (^(A),*) = .) = P (D , n = &) = P (A),r>-1 = 6) = dn-l/2 . 

This completes the proof of Prop. 16.11 

Conclusion. The following question remains: does there exist a positive-rates PCA 
which is non-ergodic with a unique invariant measure? 

Let us provide some context. By definition, a PCA has positive-rates if all its prob- 
ability transitions are different from and 1 (more formally, if / : TX — > A4(S) is the 
transition function, then \/u S H v ,\/v G E, f(u)(v) G (0, 1)). It had been a long standing 
conjecture that all 1-dimensional positive-rates PCA are ergodic. In [5], Gacs disproved 
the conjecture by exhibiting a complex counter-example with several invariant measures. 
The existence of the intermediate case (unique but non-attractive invariant measure) re- 
mains open. A priori, it is not possible to perturbate model A to get a positive rates 
example. 
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